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Abstract 

The self-wetting properties of ionic crystal surfaces are studied, using NaCl(lOO) 
as a prototype case. The anomalously large contact angle measured long ago by 
Mutaftschiev [1,2] is well reproduced by realistic molecular dynamics simulations. 
Based on these results, and on independent determinations of the liquid-vapor (7lv) 
and the solid- vapor interface free energy (7sv) [3], an estimate of the solid-liquid 
interface free energy (7sl) is extracted. The solid- vapor surface free energy turns out 
to be small and similar to the liquid-vapor one, providing a direct thermodynamic 
explanation of the reduced wetting ability of the ionic melt. 
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1 Introduction 

A long known unusual property of molten alkali halide salts is their remarkable 
inability to wet their own solid surface at the melting point, giving rise, in 
particular to a partial wetting angle as large as 48° on NaCl(lOO) [1,2]. While 
partial self-wetting must generically arise when there is non-melting of the 
solid surface, as has been long known e.g. for metal surfaces [4], there are so 
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far no theoretical studies of either effects in ionic systems. In our work [3] we 
have investigated the properties of solid NaCl(lOO) near to the melting point 
in contact with a droplet of its own liquid. 



2 Method of calculation 

We carried out extensive classical molecular dynamics simulations of a nan- 
odroplet of liquid NaCl deposited on the solid NaCl surface. For NaCl, whose 
molecular model was described by the Born-Mayer-Huggins potential as pa- 
rameterized by Fumi and Tosi [5,6], the Coulomb long-range interactions were 
treated by the standard 3-dimensional Ewald method [7]. The convergence 
parameters for the Ewald summations were chosen in order to converge the 
real and reciprocal space sums to a given tolerance (1CT 5 Hartree per atom). 

We simulated NaCl bulk systems and slabs with (100) orientation, consist- 
ing of about 4000-5000 molecular units. The chosen time step was 1 fs, and 
the typical simulation time was 200 ps. After verifying that the bulk ther- 
modynamic properties of solid NaCl (melting temperature, lattice expansion, 
volume jump at melting point [3]) are very well described we simulated the 
approach of a NaCl droplet on NaCl(100) surface. 

We prepared a 7200 atom crystalline slab, 8 atomic planes thick (15 x 15 x 
4(NaCl) 4 conventional cubic cells). The slab was gradually heated up to the 
bulk melting point, keeping two atomic layers rigidly fixed at the bottom of 
the slab, however with a lattice spacing growing with T, in order to simulate a 
thermally expanding semi-infinite solid. Separately we prepared a small NaCl 
cluster consisting of 500 NaCl molecular units. The size of this cluster was large 
enough to yield after melting a near-spherical drop, and also a well-defined 
shape when the drop is deposited on the surface. The cluster was melted by 
heating above the bulk melting temperature (Tm = 1066 ± 20 K in our model 
in good agreement with the experimental value T^ p = 1074 K [3]), and then it 
was equilibrated at T=T M for 100 ps, yielding a well-defined spherical liquid 
droplet of radius 18 A. 

Finally, we placed the drop near the solid NaCl(100) slab surface, the lower- 
most atoms of the drop a distance 4 A from the surface (fig. la). Zero vertical 
velocity was assigned to the drop center of mass but the drop expanded any- 
way to touch the surface. As soon as the drop contacted the surface, attempts 
at equilibrating the system at the theoretical melting temperature (T M ) failed, 
and the whole solid NaCl slab melted very quickly, which is compatible with 
the 'fragile' nature of non-melting of solid NaCl(100) [3]. At a much lower 
temperature than Tm (1000 K) on the contrary, the substrate remained solid, 
while the drop spread slightly only at the beginning but then crystallized, 
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forming a nice "stepped" pyramid, made up of (100) facets. We succeeded 
finally in equilibrating the liquid nanodroplet on the solid surface at the inter- 
mediate temperature of 1050K, which is only slightly below T M of our model. 
This is described in the following section. 



3 Results 

The nanodroplet and the solid slab were separately equilibrated at 1050 K. 
During the first 100 ps after contact, the droplet settled down on the substrate, 
gradually approaching the final shape (fig. lb-d). In the next 130 ps, the 
droplet survived in a (not clear if metastable or unstable, but long lived) state 
without spreading appreciably (fig. le-f). At the end of our simulation, the 
droplet - substrate system looked like fig. lh. The top view of fig. 2 shows 
that the drop spread almost circularly. 

Before proceeding with further descriptions, it is important to specify the 
thermodynamics of this situation. Because we are below T M (even if slightly) 
the final equilibrium state should consist of a flat solid NaCl(100) surface, i.e. 
the nanodroplet should completely spread and recrystallize. That however will 
take a very long time. While the nanodroplet exists, it will form an external 
wetting angle #lv (fig- 1), as well as an internal angle #sl- The latter #sl is 
irrelevant here, because it depends critically on the temperature and on the 
time. In particular at short times and not too far from T M we expect 6*sl ~ 0. 
The external angle #lv is instead significant, as it should equal the macroscopic 
wetting angle measured in the bubble experiment [1,2]. 

Assuming quasi equilibrium for the liquid nanodroplet in the sense specified 
above, this angle will obey Young's equation (fig. 3) 

7sl + 7lv cos 6> lv = 7sv (1) 
where the 7's are the interface free energies. 

To determine the external wetting angle 6*lv of the nanodroplet we analyzed 
100 configurations in last 100 ps. The instantaneous atomic positions were 
plotted in cylindrical coordinates (r and z, where r is parallel to the surface), 
and from the profile of the drop, we determined the best approximation to a 
portion of a sphere, by determining the center position and the radius. The 
contact angle follows immediately by simple geometry from these two quanti- 
ties. Our best estimate slightly below Tm is 6?lv = 50°± 5° which is in excellent 
agreement with the experiment value at the melting point (48°) [1]. At the end 
of the simulation the internal solid-liquid interface was still relatively sharp 
and flat, consistent with our assumption 9 SL ~ 0. 
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The thermodynamical reason leading to the large observed #lv can now be 
analyzed. As eq. (1) shows, a large value of #lv may arise either due to a large 
value 7sl, or to a small 7sv, relative to 7lv or to both. These free energies at 
T m T M are all presently unknown. While simulations are presently under way 
to calculate 7sv [3], we can already anticipate on purely physical arguments 
that both factors namely, a large 7sl and a small 7sv, will indeed occur on 
NaCl(lOO). The liquid and the solid differ enormously, in density and other 
properties, and this does suggest that 7sl has no reason to be small, as it was 
instead in the case of metals. The solid surface free energy is instead small 
due to a good cancellation of the Coulomb potential outside the (100) surface, 
and also to a good vibrational free energy of the hot solid surface [3]. 



4 Discussion 



We can now try to relate further the wetting angle #lv to more microscopic 
quantities. Di Tolla et al. [4] showed that for a non-melting metal surface, the 
angle #lv could be related through a simple model to the surface spinodal 
temperature Ts, that is the maximum temperature above Tm to which the 
non-melting surface can be overheated. As discussed in a companion paper to 
this [3], the very same model must be improved in order to describe NaCl, 
owing to the larger differences between solid and liquid. The free energy gain 
upon crystallizing the surface at T = T M , A7oo = (7 LV + 7sl) ~ 7sv is related 
to the partial contact angle through Young's equation (1) 

cos# LV = l-^ (2) 
7lv 



If we plug in our calculated value of 7lv — 104erg/cm 2 and = 50°± 5°, we 
obtain A 7,^ = 37 erg/cm 2 . The large value of A7oo represents the "binding 
energy" of the solid-liquid and liquid-vapor interfaces, to form a single one, 
namely solid-vapor. Assuming a standard form [8] 

1^(1) = -^ cos ^e"V€ (3) 

2 CL 



for the interface interactions, it leads to a connection [3] between A7oo and 
the solid surface spinodal temperature Ts, which is the temperature where 
the solid surface becomes mechanically unstable, in the form: 

A 7oo ~pLa(Ts/T M -l), (4) 
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With our value of A7oo, and L = 4.813 x 10 9 erg/gr, a = 5.9 A we predict 
T s = 1210 K, very close to that seen in simulations. [3] 



5 Conclusion 

In conclusion, molecular dynamics simulations indicate that metastable liquid 
NaCl nanodroplets on solid NaCl(lOO) can exist for some time close to Tm- 
They are found to exhibit an unusually large partial wetting angle, whose 
value is in agreement with the large macroscopic wetting angle observed in 
the bubble experiment. Thermodynamically this appears to be caused by an 
unusual small value of the solid-vapor interface free energy, as well as by a 
large solid-liquid interface free energy. The free energy differences extracted 
via Young's equation can also be connected well with other properties, such 
as the solid surface spinodal temperature. 
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Fig. 1. Time evolution of the liquid NaCl drop on NaCl(lOO). Dark and light circles 
stand for the Na + and CI - ions respectively. 
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Fig. 2. Top view of the liquid droplet on the NaCl(lOO) after 230 ps. Dark and light 
circles stand for the Na + and Cl~ ions respectively. 




Fig. 3. Sketch of the liquid drop partially wetting a solid substrate, showing the 
balance of the forces acting at the interfaces. 
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